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Abstract. Frequency separations used to infer global properties of stars through 
asteroseismology can change depending on the strength and at what epoch of the stellar cycle 
the p-mode frequencies are measured. In the Sun these variations have been seen, even though 
the Sun is a low-activity star. In this paper, we discuss these variations and their impact on 
the determination of the stellar parameters (radius, mass and age) for the Sun. Using the data 
from maximum and minimum activity, we fitted an age for the Sun that differs on average by 
0.2 Gyr: slightly older during minimum activity. The fitted radius is also lower by about 0.5% 
for the solar effective temperature during minimum. 



1. Introduction 

The p-mode oscillation frequencies vary with the activity cycle of the Sun pQ and other solar- 
like stars [2]. However, as observed in the Sun, these temporal variations present differences 
in amplitude, phase, and frequency among individual modes. Thus the frequency separations 
computed will be observed to vary with time. Comparing these frequency separations with those 
from stellar models will lead to differences in the fitted global parameters (radius, mass, age) of 
the star. Generally, when the frequencies are measured we do not know the phase of the activity 
cycle, unless the star can be measured continuously for long periods of time. We should then 
take into account when inferring the stellar parameters, that these values may depend on when 
the data was obtained. 



2. Analysis of time-series observations and determination of frequencies 

We analyzed 5202 days of velocity observations collected by the space-based instrument GOLF 
onboard SoHO spacecraft, covering more than 14 years between 1996 and 2010 with an overall 
duty cycle of 95.4% [3]. This dataset was split into non-independent, contiguous 365-day 
subseries with 91.25-day overlap and their associated power spectra fitted to extract the mode 
parameters [4] using a standard likelihood maximization function (power spectrum with a x 2 
with 2 d.o.f. statistic). The formal uncertainties in each parameter were then derived from the 
inverse Hessian matrix. 

The individual frequency separations measured with 365-day time series were obtained with a 
precision of about 0.07 /iHz. Figure [1] (left panel) shows the temporal variations of the individual 



I = and I = 2 mode frequencies averaged between 2200 and 3300 /^Hz (the reference values 
being taken as the average over 1996-1997). It is clear that the individual Sun-as-a-star p- 
mode frequencies have different temporal variations pQ, which are consistent between radial 
velocity (GOLF) and intensity VIRGO observations [5]. The averaged large (Az^) and small 
(5i>) frequency separations also show significant temporal variations with solar activity. For 
example, Figure [1] (right panel) shows that the small frequency separation 5vq2 varies from 
peak-to-peak by about 0.2 ± 0.02^Hz over the solar cycle, which is presumably consistent with 
being signatures from surface effects. For a broader perspective, it is very important to note 
that the Sun is considered a low-activity star, and many solar- type stars could exhibit much 
larger variations with activity cycle. 

3. Solar global parameters 

3.1. Fitting strategy 

In order to determine the global model parameters of the Sun P, we match the observational 
data O to the observables M from stellar models [HE] which are characterized by P. P comprises 
the mass M, age A, initial hydrogen (or helium) mass fraction Xq, initial heavy element mass 
fraction Zq, and the mixing- length parameter a (from the standard mixing-length theory of [8]). 
The parameters that best describe the data are obtained by minimizing a \ 2 function: 



where there are i = 1,2, ...,M independent observations and a is the observational error. We 
use the Levenberg-Marquardt algorithm (LM) to minimize Eq. Q] [9[ [TUl [TTT [TZj [T3] . 

The global objective of this work is to investigate the effect of the shifted frequency values 
on the determination of P for any star with solar-like oscillations. So we must define a coherent 
and consistent method to fit the observational data to obtain the set of best-fitting parameters 
not only for the Sun, but for other stars too. Our strategy is to first fit the non-seismic and the 
average seismic quantities to obtain P to a reasonable range, and then using P as initial guesses, 
proceed to use the individual frequency separations (large and small) to obtain the final set of 
best-fitting parameters. 

LM is a local minimization method and can therefore be sensitive to initial conditions. To 
ensure that we obtain a global minimum we search for the best parameters by initializing the 
minimization with different parameter values. We set the initial M as 1.0 Mq, and the initial 
A varies between 4 and 6 Gyr in steps of 0.5 Gyr. Additionally, we hold the parameter Xq fixed 
during each minimization, and instead repeat the process for three values of Xq: 0.69, 0.71, 
0.73. This gives in total 15 minimizations using the input data at maximum activity and again 
at minimum activity. 



From the 15 best-fitted P, we selected those sets where x 2 < 29.1 (4 fitted parameters, 19 data 
points to give 14 degrees of freedom at the 1% confidence level) and where the fitted T e g fell to 
within la of 5777 K (T e g was not a constraint in the second part of the fitting strategy). These 
sets are represented in Fig. [21 where we show the values of age versus mass on the left panel, and 
radius versus effective temperature on the right. The dark/lighter filled circles are the results 
for fitting during maximum/minimum activity. 

Fig. [2] left panel shows an offset in the fitted P between minimum and maximum activity. To 
quantify this offset, we fitted a linear function to each group of results. The difference between 
the fitted ages at opposite activity phases is 0.2 Gyr. Alternatively if we take 4.8 Gyr as the 
solar age, the mass is fitted as a 1.0 Mq star at minimum, and 1.03 Mq at maximum. 
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3.2. Results 
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Figure 1. Left: Frequency shifts (/uHz) of the I = (• ) and 1 = 2 (■) modes measured from 
GOLF data. Right: Temporal variations of the averaged small frequency separation <5z/Q2 (/xHz). 



The right panel of Fig. [2] shows the fitted R and T e g for these same models (neither were 
input constraints). We find that at minimum activity, the Sun's T e g is about 50 K cooler for a 
given R, or about 0.5% smaller in R for a given T e g. 

4. Conclusions 

Measuring the oscillation frequencies during different phases of a stellar activity cycle will lead to 
different values of the individual frequencies. Thus, estimates of the large and small frequency 
separations will be different depending on the observing period. Using 365-day subseries of 
GOLF data, we determined the individual frequency separations to about 0.07 //Hz. 

We used stellar models to determine the best-fitting parameters using the observations taken 
at minimum and maximum activity. We found that the age of the star is on average 0.2 Gyr 
older using the values from minimum activity, or for a given age, the mass is about 2% larger 
using the data at maximum activity. We also found a small discrepancy in the fitted radius 
and effective temperature. At minimum activity, the Sun is about 50 K cooler for a given fitted 
radius, or the radius is about 0.5% smaller for a given effective temperature. 

Although we still need to study whether these differences in fitted values are detectable given 
the expected uncertainties, we note that [13] quote an error in the age of the planet-hosting star 
HAT-P-7 of 0.26 Gyr, of the order of this detected change. The Sun, however, is considered a 
low-activity star, and much stronger activity cycles on other solar-type stars have already been 
detected [2]. With long time-series data such as those from CoRoT or Kepler, we will not only 
detect stellar activity cycles, but we will obtain very high precision on the seismic data, that the 
uncertainties in the fitted parameters will be smaller than the fitted changes in these parameters. 
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Figure 

separations calculated at maximum (black) and minimum (grey) activity, 
corresponds to a minimization. The lines are linear fits to the results 
radius and effective temperature for the same minimizations. 
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